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A simple model describing depolarization channels with zero-bandwidth environment is presented 
and exactly solved. The environment is modelled by Lorentzian, telegraphic and Gaussian zero- 
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bandwidth noises. Such channels can go beyond the standard Markov dynamics and therefore can 
illustrate the influence of memory effects of the noisy communication channel on the transmitted 
information. To quantify the disturbance of quantum states the entanglement fidelity between 
arbitrary input and output states is investigated. 
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I. INTRODUCTION 



One of the most important features concerning quantum communication is the capacity 



or the fidelity of quantum information transmitted in noisy quantum channels The key 
factor limiting the possibilities of communication using quantum states, is an environment- 
induced noise. Uncontrolled interaction between the environment (E), and the transmitted 
quantum state can essentially affect the state and in consequence lower the communication 
capacity of the information channel 

aa. 

For qubits, a well known class of quantum noisy channels consists of depolarizing channels 
P). Input information of such channels is stored in a density operator g in . Such channels 
can be characterized by a probability p that the quantum information is distorted, and with 
a probability 1 — p that the information remains intact. In the simplest case of a single 
qubit transmitted through the noisy channel the influence of noise is usually decomposed 
into three interaction channels. Bit error channel a x flipping the values of bits: |0) i— > |1), 
|1) i — > |0); phase error channel b z flipping the phase: |0) i— > |0), |1) i— > — 11); and phase and 
bit error channel a y flipping both: |0) i— > |1) i— > — i\0). 

The influence of these interaction channels can be written as an incoherent combination 
of three unbiased terms generating bit flip errors and phase flip errors in the form given by 

Q: 



gout = (1 -PjQin 
p 

This depolarizing channel is just an example of a general quantum channel characterized 
by a trace-preserving, linear map $ : £> in i— > g ont . The influence of the environment on 
the quantum state g in can be represented in terms of a Kraus decomposition: $(ft n ) = 

Ertftn^J, where E r ^Jt = 1 fl- 

The unbiased depolarizing channel given by Eq. (JTJ) , corresponds to Kq = y/1 — p 1 , 

It is not difficult to derive the expression corresponding to Eq. ([!}, from a unitary evolu- 
tion involving an extended Hilbert space TC®TCe, of the input qubit and an additional qubit 
(E), consisting the environment degree of freedom. Although mathematically correct, such 
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a derivation has no simple or direct physical realization in terms of a realistic noise. In fact 
in most known cases, the environment is much more complex, and the additional fact that 
the three interaction channels corresponding to bit error, flip error and phase error do not 
commute, leaves the incoherent addition of these channels in the formula (Q) questionable. 

A more realistic approach to depolarization channels with bit errors requires a better 
understanding of the physics involved in the system-environment interactions. We shall 
assume that the system-environment interaction is characterized by a model Hamiltonian 
H. A unitary evolution of this combined system leads to a Kraus decomposition described by 
a time dependent map $t : g- m i— > g(t). This means that one should have a time- dependent 
p(t) such that at the input time t = 0, p(0) = 0. Because of this the general property of the 
map: &t\t=o — 1? expresses the continuity at the origin and hence for all time. 

In general the problem of finding the evolution of the state interacting with the environ- 
ment is very difficult, therefore the environment- induced noise is modelled using various sim- 
plified approaches including several assumptions such as the Markov property. In such case 
it is assumed that the evolution of the state at a given instant is fully determined by the state 
at that instant, so the process has no "memory" of its past. The Markov property means 
the quantum channel is such that for an infinitesimal time interval: p(At) ~ At. In this 
case the channel map generates a completely positive dynamical semigroup: $t o $ t / = 
for t,t' > 0, which defines a Markovian dynamics. 

As it is well known, for such Markovian maps we can transform the time-dependent Kraus 
decomposition (JTJ) into a local Lindblad equation pj 

^ = &<«). P) 

with the initial condition g(0) = g- m , and where the Lindblad superoperator £ can be derived 
from the Kraus operators. 

In this paper we present a simple model describing an evolution of a quantum state 
interacting with the environment and study various properties of the affected state. A 
physical picture behind the algebra can be for example a randomly fluctuating magnetic 
field acting on an electron's magnetic moment, or a thermally fluctuating birefringence of a 
single mode fiber transmitting a polarization state of a single photon. 

Using our model we will justify the assumption and determine the conditions for 
its validity. We will show that the disturbed output state has the form (JTJ) in the infinite 



interaction time limit only if the disturbance is unbiased and acts identically in all bases. In 
the general case the evolution leads to a different final state. It will also turn out that the 
simple model leads to a nontrivial evolution not obeying the Markov property. Therefore, 
our simple approach will lead us out of the no-memory approximation regime. Our model 
shows that in many cases one cannot fairly neglect the memory effects of the interaction 
and consequently the Markov property is not always valid 0,0]. 

The paper is organized as follows: in Sec. |H]we present the model of a qubit in a presence 
of an external noise. In Sec. IIHI we show that if the environment of the channel consists 
of zero-bandwidth noises one can calculate exactly expressions determining the evolution 
of a single qubit in such channel. Sec. I1VI contains several examples of zero-bandwidth 
noises that can affect the input state. We show that in general the decoherence channel can 
not be regarded as an incoherent superposition of independent interactions. Special cases 
involving Markovian and exactly soluble non-Markovian dynamics are derived. In SecfVlthe 
efficacy of the quantum channel is quantified by the fidelity between the output state and 
the input state. An appropriate measure for assessing the fidelity of a mixed input state is 
the entanglement fidelity |lfj , which is the maximum fidelity of states being purifications of 
the input mixed state p and the output state $(p). Finally Sec. EU concludes the paper. 



II. INTERACTION MODEL 



The interaction of a qubit with an environment inducing random bit errors will be de- 
scribed by the following Hamiltonian 

H — r ■ & — xa x + ya y + za z , (3) 

where the three components r« = (x,y,z) will be uncontrolled "noisy" parameters charac- 
terizing the fluctuations of the environment. 

We will assume that are independent random variables. This situation is a simplifica- 
tion of a general very difficult problem with rj-(i) being time-dependent stochastic processes 
with arbitrary autocorrelations: (ri(t)rj(t')) = 5^ Ai(t,t'). In these applications, the au- 
tocorrelation functions Aj(t, £'), have usually a Fourier limited spectrum with an effective 
bandwidth 7 characterizing the environment noise. Even for the simplest form of the auto- 
correlations, the exact solution of the full time-dependent problem involving more that one 
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Ti(t) noise is not known. However our simple noise model can illustrate several properties 
of various completely positive maps. 

In the proposed scenario, the evolution of a single qubit given by the von Neumann 
equation has the following form 

~jf = i Ti t^' S] = A>£ (4) 

i 

The appearing Liouville superoperator Cq describes a unitary rotation of a qubit defined 
by the coefficients rj. The solution involves a stochastic averaging with respect to the 
environment. As a result it can be written in the following compact and formal form: 

e(t) = {T<£ ia£oW )Q**- (5) 

There is no useful formula that can handle the chronological time-ordering of three Pauli 
matrices, and allows an exact stochastic averaging over the environment noises. There are 
however special cases when the exact solution of Eq. (J3J) can be obtained. We know of 
three cases. Case one involves an arbitrary stochastic noise and only one rj. In this case 
the chronological ordering plays no role. In case two, fluctuations are Gaussian and the 
bandwidth characterizing the environment noise is infinite 7 = 00 (white noise). In this case 
an exact average of Eq. (j5J) exists, and a Lindblad equation (J2J for the channel map can 
be derived. Case three, the one investigated in this paper, corresponds to arbitrary random 
fluctuations of r, with the environment described by a zero-bandwidth environment noise: 
7 = 0. In this case all autocorrelations Aj(t, t'), become time-independent, the chronological 
product plays no role, and an exact averaging over the environment with arbitrary statistics 
can be performed. 

III. EXACT SOLUTION WITH ZERO-BANDWIDTH 

In order to find the evolution of a state under a time-independent Liouvillian C , we first 
find its eigenstates: 

C A ■ & = iriAjl&i, frj] = -2riAje ijk a k = AA ■ <r, (6) 
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hence we obtain a set of equations: 



-IriAjtijk = \A k . (7) 

The solutions exist only for a set of eigenvalues A G {0, 2ir, —2ir}, where r = \J x 2 + y 2 + z 2 . 
For this set of eigenvalues we find the corresponding eigenvectors: for A = we have A = r, 
and for A± = ±2ir we have A± = (=pyr — xz, ±ixr — yz, x 2 + y 2 ). 
In order to determine the evolution of an arbitrary initial state 

g iQ = l(i + a-d-), (8) 

it is helpful to decompose it into the calculated eigenvectors. Pauli operators &i written in 
the calculated eigenbasis have the following form: 



x xz(A + + A_) - iyr(A + - A_) \ „ 
r 2 2r 2 (x 2 + y 2 ) J 

f y_ yz(A + + A_) + ixr(A+ - A_) ^ „ 
°v ~ y r 2^ 2r 2 (x 2 + y 2 ) J ' * 

At this point it is easy to find the action of the evolution operator exp(£ot) on the given 
input state. We simply multiply the eigenvectors appearing in our decomposition by the 
proper factors exp(At) with corresponding eigenvalues A. This yields: 



e fin — 2 



. x xz(e 2irt A+ + e" 2ir *A_) - iyr(e 2irt A + - e~ 2irt A_)\ 
° fa ^f^ Ao 2r 2 (x 2 + y 2 ) J" 7 



y yz(e 2trt A + + e~ 2irt A_) + ixr(e 2 ' trt A + - e~ 2irt A_) 



+a y ~^A - 



2r 2 (x 2 + y 2 ) 



+n ,(± Ao+ ^ A % + r r A - )^]. m 

The above formula expresses the state of the qubit evolving under the action of the 
Liouvillian defined by the arbitrary vector r. 
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In our model of the noisy channel, the interaction between the environment and the qubit 
can be described by a randomly chosen vectors r. Therefore, to model the evolution of the 
qubit under the influence of the environment-induced noise we will average the obtained 
output state over all possible realizations of the dynamics characterized by arbitrary vectors 
r. For simplicity we will be interested in an averaged evolution of the qubit with an even 
in r probability distribution p(r) = p(—r). In this case a non- vanishing contribution to the 
averaged output state will come only from the symmetric part of the expression (fTU|) : 



e Cot g in 



sym 



y2 _|_ ^2 

1 + a x a x ( — H — cos 2rt 



<£ 2 _|_ £ 2 

+ a v°y [ ^2 ^ ~2 — cos 2rt 



2^ -\- y 2 
+ a z a z | — H - — cos2rt 



(11) 



From the above formula it follows that the initial state (JHJ) evolves into an averaged output 
state 



£?out(£) = - (i + aihi{t)<Ti) . 



(12) 



The dynamics of the output state is completely described at any time by a set of time- 
dependent functions Aj(t), which have the following form: 



Ai(t) = 1 - 2 j d 6 rp(r) - -|J sin 2 rt. 
The output state can be also equivalently represented in terms of Kraus operators K r : 



(13) 



where 



eout(t) =J2^r(t)Q in kl(t), 



(14) 
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Ko = ^1^1 + A x + A y + A z = tJ J d 3 rp(r)cos 2 rt, 



K\ = ^xV 1 + A x -A y -A z = a x J J d 3 rp(r)^-sin 2 rt, 



K 2 = ^a y y/l - A x + A y - A z = a y J J d 3 r p(r)^ sin 2 rt, 



K 3 = ^a z ^/l-A x -A y + A z = a z] j j d 3 rp(r)^sin 2 rt. (15) 

These exact expressions for the Kraus operators, are the main result of our investigations. 
Before we discuss various statistical models, we note that an unbiased incoherent addition 
of bit-error channels in most cases is not justified. The formula above shows that the 
time-dependent A«(t) functions couple in a highly nontrivial way the three channels. The 
simplified expression (0), does not reflect this complicated entanglement between various 
bit-error channels. 



IV. EXAMPLES OF A NOISE 
A. Markov noise 

Consider a simple case of a completely positive map determined by a Lorentzian proba- 
bility distribution 

p(r > = h (f^wm 

(16) 

characterized by a width T. The corresponding Kraus operators (fT5|) in this case read 



S 



#0=1 

Ki = a x 



K 2 = <3v 



# = a. 



1 + exp(- 


-rt) 


2 


/l — exp( 


-n) 


2 


/l -exp(- 


-n) 


2 


/l — exp(- 


-rt) 



(17) 



The resulting dynamics of such a channel is: 



1 + exp(-rt) 1 - exp(-rt) 
Q{t) = g + 6 



X [(^(O)^ + fryQ^fry + <T Z ^(0)(J 2 



(18) 



Let us note that this expression is equivalent to Eq. (JTJ if the probability is time-dependent 
i.e., 



p(t) 



1 - exp(-rt) 



In the steady state p(oo) = |, the quantum channel reduces to a very simple expression 



(19) 
ion j^: 



1 



1 



1 



1 



Pout = 2 I Pin + ^VxQinVx + ^yQm^y + ^zQin^z 



(20) 



One can easily check, that the infinitesimal time evolution of the last three Kraus opera- 
tors is: Ki(At) ~ \/At<jj. This behavior is typical for a diffusion process and consequently 
the state evolution clearly obeys the Markov no-memory property, and as a consequence has 
the form of the Lindblad equation (J2J): 



dg 
dt 



Cg{t) 

r 

~2 



Q - g {^xQcr x + CTyQCTy + a z ga Z/ 



(21) 



B. Telegraphic non-Markov noise 

Now, let us assume, that the noise introduced to the system is a random telegraphic noise 
[ll| . so that the disturbance of the qubit induced by the environment is discrete, and jumps 



between two values ±a. For concreteness we consider the following probability distribution: 



p(r) = ~ [S(x - a) + <5(x + a)] S(y)5(z). (22) 



In this case the Kraus operators ([15]) equal: 



Kn = iV cos 2 at, 



Ki = d x V sin 2 at, 

= £ 3 = (23) 
and the disturbed qubit at instant t is in the state: 

g(t) = cos 2 at g(0) + sin 2 at a x g(0)a x . (24) 

The periodic result is very straightforward, however it reveals something interesting. Al- 
though our model is quite simple, it leads to non-trivial dynamics, which becomes apparent 
when we analyze the evolution of the density operator One can easily find that the 
time-evolution is given by a non-local in time Lindblad equation: 

<W) " 2 ft ds[g(s)-a x g(s)a x }. (25) 



dt 2 J( , 

This shows that the time dynamics of g(t) at the given instant t depends not only on 
the state at this instant, but also on the state at all earlier times. This behavior can be 
seen already from the form of the Kraus operator, which for the infinitesimal time evolves 
as: Ki(At) ~ Ata x and such an evolution characterizes non-Markov processes with zero- 
bandwidth 0| . 

From this example we conclude that our model in general does not obey the "no-memory" 
approximation and the evolution of the state is non-Markovian. One could also think of 
studying multi-dimensional telegraphic noise, however in this case the analysis becomes 
much more complicated and there is no simple, linear integral kernel as the one in Eq. (|25jl. 

C. Gaussian noise 

Although in general the expressions ()15|) are not analytically integrable, one can find 
explicitly the Kraus operators in the asymptotic steady-state limit t — ► oo. In this limit 
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FIG. 1: Kraus operator's coefficients ki (denned via the relation: Ki = ki&i) for a Gaussian 
probability distribution p(r) = — ===== exp (— §r — |r — ) , where on the upper figure d x = 
d y = d z = 1 and on the lower = 1, d y = 2, = 3. 

the square of rapidly oscillating trigonometric functions appearing in the integrals (jl5)l can 
be approximated by their average value ~. In the simplest case of an arbitrary, spherically 
symmetric probability distribution p(r) the coefficients ki of the Kraus operators (defined 
via the relation K~i = ki&i) equal k = k\ = k 2 = ks = -j= and consequently the output 
quantum state reads as in Eq. ([20)1 . 

This result reproduces the steady-state Markov limit justifying its validity. However it 
is valid only when the probability distribution p(r) is spherically symmetric i.e., the three 
channels are unbiased. In general the input state evolves to a different limit. 

In the Figure ^ we have shown the numerically calculated evolution of the Kraus opera- 
tor's coefficients for a Gaussian probability distribution 

p(r) =7p|n exp (-S-5ri) (26) 

for the following two cases: when the probability p(r) is spherically symmetric with d x = 
dy — dz — 1 (upper plot) and for the asymmetric distribution with d x — 1, d y — 2, 
d z — 3 (lower plot). It is seen that after some characteristic time, the state becomes 
stationary, however the limit depends on the characteristics of the probability distribution 
p(r). The example discussed above provides an illustration of a unbiased and biased Gaussian 
depolarization channels. 



11 



V. FIDELITY FOR MIXED INPUT STATES 



In order to judge the quality of a communication channel and the role of the introduced 
noise one needs a tool to investigate the state disturbance during the transmission. To 
quantify the influence of the external noise onto the transmitted quantum state we use an 
entanglement fidelity measure defined as the following overlap between the input and output 
density matrix J|: 

HQ^ Qont) = (Trjv^outftt}) • (27) 

This fidelity is in general very difficult or impossible to calculate. For an arbitrary input state 
of a single qubit given by Eq. (jHJ, and with an arbitrary spherically symmetric probability 
distribution p(r) characterizing the external noise this fidelity can be calculated exactly and 
is equal to: 

HQ^ Qont) = l(t+ Vx(l~a 2 )) , (28) 

where £ = 1 + a^.K x + a 2 A y + a 2 z K z and x — 1 — a l^-l ~ a y-^y ~ a z-^x- For pure states (a = 1) 
the formula simplifies to T = |. On the other hand it is not very surprising, that the 
maximally mixed state (a = 0) remains unchanged under the influence of the noise, while 
the communication fidelity decreases with increasing purity of the input state. 

Using the same approach one may study also the evolution of multidimensional systems. 
Of course the general expressions become very complicated, even when we consider a two- 
qubit Hilbert space, however it is possible to find some compact solutions for special cases 
of pure states. 

Consider an arbitrary two-qubit initial pure state: 

l^n) = a\ ID + b\ !<->) + c\ <->J) + d\ <-«->>. (29) 

Using the same approach as above, one can calculate that for the independent disturbance of 
each mode with the same type of noise characterized by the spherically symmetric probability 
distribution p(r) the fidelity of the transformation is: 
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FIG. 2: Fidelity measure as a function of time t for the input state Q[ n — 
\ (I tl)( tt I + I I) + y (I tt)( | + | tt I) for several values of the parameter m. 

From bottom to the top, respectively: m = 1, m = 0.9, m = 0.7, m = 0.4, m = 0. 



•F(|*>in, &ut) = (^{Vl^inX^inleoutl^inX^ml} 
(^m|^out|^in) 



( 1 ^) 2 -4A^|6c-ad| 2 . (30) 



What is interesting in the above result is that the fidelity is the highest for separable states 
(for example for a = b = and any c, d) and it drops down when the input state becomes 
more entangled. 

Another compact result can be found for the following two-mode mixed input state: 

ftn = ^(ItDUII + l— ><<-><-> I) 
777 

+y(itt)< — I + !<-><-> XXII)- (31) 
With a similar analysis one obtains the fidelity measure given by: 
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F (g- m , q, 



but; 



1 + A 2 Z + m 2 (A 2 X + A 2 y )+VT 



x^/(l+A2) 2 -m2 (A2+Ag) 



(32) 



In the Figure El we have plotted the dynamics of fidelity for several parameters m and the 

3 2 

unbiased Gaussian probability distribution p(r) = 7r~2e~ r . We find a not very surprising 
result, that the transformation fidelity is a decreasing function of the purity of the input 
state. For m = (which of course does not yet correspond to the maximally mixed state) 
the fidelity is the highest, while for the pure state [m = 1) the fidelity is the lowest. 



VI. CONCLUSIONS 



Uncontrolled interaction between the environment and the transmitted quantum state 
can essentially affect the state and in consequence lower the communication capacity of an 
information channel. Several ideas has been put forth to overcome this problem. One of the 
m08t proving is the ase of so-cal.ed deeoherence-free Spaces Ql3. T hlS ,dea eaa be 
applied when the noise present in the system is correlated between consecutive uses of the 
communication channel Q|, however this is not always possible and therefore one needs a 
careful study of the properties of various types of noise and their influence on the quantum 
state. 

In this paper we have introduced a dynamical model of interaction between the quantum 
state and its environment and shown that although based on simple assumptions, it leads to 
non-trivial solutions. Using the model we have analyzed properties of zero-bandwidth noise 
with Lorentzian, telegraphic and Gaussian distributions and have shown that only the first 
of them obeys the Markov property, while the others exhibit memory effects and are non- 
Markovian. Our approach allowed us to solve a simplified version of a general problem when 
the noise is an arbitrary time-dependent stochastic process, whose solution is not known. 
We have calculated transformation fidelities for a collection of input states and analyzed 
their dynamics according to their entanglement or purity. 
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